Resonant scattering in graphene with a gate-defined chaotic quantum dot 



Martin Schneider and Piet W. Brouwer 
Dahlem Center for Complex Quantum Systems and Institut fiir theoretische Physik, 
Freie Universitdt Berlin, Arnimallee I4, 14195 Berlin, Germany 

We investigate the conductance of an undoped graphene sheet with two metallic contacts and 
an electrostatically gated island (quantum dot) between the contacts. Our analysis is based on the 
Matrix Green Function formalism, which was recently adapted to graphene by Titov et al. [Phys. 
Rev. Lett. 104, 076802 (2010)]. We find pronounced differences between the case of a stadium- 
shaped dot (which has chaotic classical dynamics) and a disk-shaped dot (which has integrable 
classical dynamics) in the limit that the dot size is small in comparison to the distance between 
the contacts. In particular, for the stadium-shaped dot the two-terminal conductance shows Fano 
resonances as a function of the gate voltage, which cross-over to Breit-Wigner resonances only in 
the limit of completely separated resonances, whereas for a disk-shaped dot sharp Breit-Wigner 
resonances resulting from higher angular momentum remain present throughout. 
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I. INTRODUCTION 

With its remarkable electronic properties, graphene, a 
single layer of carbon atoms as they occur in graphite, 
continues to be the subject of theoretical and experi- 
mental research.—^— Graphene owes its special electronic 
properties to a Dirac-like dispersion in the absence of 
impurities, with the Fermi level at the Dirac point in 
the absence of external doping. Although the den- 
sity of states vanishes at the Dirac point, transport in 
clean undoped graphene is possible through evanescent 
modes, giving rise to a finite conductivity Ae'^/nh at zero 
temperature4i^— (At finite temperature, the conductivity 
of clean undoped graphene is formally divergent if effects 
of electron-electron interactions are neglectedjii"— ) 

If the Fermi level ep is moved away from the Dirac 
point e = 0, e.g., by application of a gate voltage, a finite 
concentration of charge carriers is created, for positive as 
well as negative ep. It is for this reason that generically 
charge carriers in graphene can not be confined electro- 
statically, or that the addition of charged impurities to 
otherwise undoped graphene leads to an enhancement of 
the conductivity, rather than a suppression. (For this 
effect to occur it is essential that the charged impurities 
are at a finite distance from the graphene sheet, so that 
they do not couple the Fermi pointS;- ^^ ) 

Recently, Bardarson, Titov and one of the authors 
showed that there is an exception to the rule that elec- 
trons in graphene can not be confined with the help of 
gate voltages, if the shape of the gate-defined region is 
such, that the Dirac equation is separable.^" An example 
of a gate-defined region (or "quantum dot") for which 
the Dirac equation is separable is a disk-shaped dot with 
a uniform potential inside the disk. In this case, elec- 
tronic states have a well-defined orbital angular momen- 
tum. Electrons with the minimal angular momentum are 
not confined, whereas electrons with a higher angular mo- 
mentum are bound to the quantum dot. The confinement 
is perfect if the graphene sheet surrounding the dot is un- 
doped (with £p =0); otherwise the bound states become 



resonances If the shape of the quantum dot is such, 
that the Dirac equation is not separable, as is the case for 
a stadium-shaped gated region, no bound states exist. 

In a semiclassical picture, the difference between in- 
tegrable and non-integrable shapes can be understood 
with the help of Klein tunnelingj^i^i An electron inci- 
dent on an electrostatic barrier in graphene is transmit- 
ted perfectly at perpendicular incidence, while the reflec- 
tion probability sharply increases away from perpendic- 
ular incidencei ^°'^^i^^ For a quantum dot with a shape, 
such that the classical dynamics is integrable, the angle of 
incidence is fixed, or only a limited set of angles appear in 
an electron's orbit, which allows for the confinement j^^i^^ 
On the other hand, with non-integrable classical dynam- 
ics, every electron eventually will approach the interface 
at perpendicular incidence and escape. (While this semi- 
classical argument explains the qualitative difference be- 
tween integrable and non-integrable gate-defined quan- 
tum dots, it fails to account for details of the bound or 
quasi-bound states, because it does not take into account 
the Berry phase, which leads to half-integral angular mo- 
menta in graphenei^) 

As shown in Ref. '20I, the existence of bound states can 
be revealed through a two-terminal measurement of the 
conductance of an undoped graphene sheet with a gate- 
defined quantum dot (see Fig. [T|), which shows narrow 
resonances for the bound states, and wider resonances for 
the quasibound states. The difference between integrable 
quantum dots (with bound states) and non-integrable 
quantum dots (without bound states) becomes most pro- 
nounced in the limit that the distance L between the two 
metallic contacts becomes large. However, a detailed in- 
vestigation of this limit was not possible in Ref. [13 be- 
cause of limitations of the numerical approach required 
to study the non-integrable case. 

Recently, the problem of gate-defined quantum 
dots in graphene was revisited by Titov, Ostro- 
vsky, Gornyi, Schuessler, and Mirlin.'^" These authors 
adapted the matrix Green function method, origi- 
nally developed by Nazarov in the context of meso- 
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FIG. 1: (Color online) Stadium quantum dot (parameters R, 
a) in a rectangular graphene sample of dimensions L x W. 
The sample is attached to metallic leads on the left and right 
side, defining a two-terminal geometry. 



scopic superconductivity^ to two-terminal transport in 
graphene. One of their main results is a relation be- 
tween the two-terminal conductance of a rectangular 
graphene sheet containing an arbitrary "scatterer" and 
the scatterer's T-matrix. A gate-defined quantum dot is 
a special case of such a scatterer. Using this method, 
Titov et al. were able to give an analytic expression 
for the conductance of an otherwise undoped graphene 
sheet with a disk-shaped quantum dot, reproducing the 
numerically obtained conductance of Bardarson et al^ 
The usefulness of the matrix Green function method has 
also been demonstrated by its application to disordered 
graphene j^i^ 

The purpose of the present article is to find out what 
information the matrix Green function method makes 
available if we consider undoped graphene with a non- 
separable quantum dot. This requires a numerical cal- 
culation of the T matrix, since no analytical results are 
available for a quantum dot for which the Dirac equation 
is not separable. We here show that the numerical eval- 
uation of the T matrix can be done very efficiently and 
that the remaining steps in the matrix Green function 
formalism can still be carried out analytically. In partic- 
ular, once the T matrix of the quantum dot is known, 
the L dependence of the conductance resonances can be 
found without additional numerical effort. This allows 
us to obtain information about the constitution of the 
resonances deep into the regime in which resonances are 
separated, which could not be addressed with the generic 
numerical algorithm used in Ref. [l^. We apply the for- 
malism to a stadium-shaped quantum dot. Our main 
result is that the resonance lineshapes are described by 
the Fano resonance formula in the case of the stadium- 
shaped {i.e., non-separable) quantum dot, while they are 
Breit-Wigner resonances for the disk-shaped dot. 

This paper is structured as follows: In Section|lI]we de- 
scribe the matrix Green function method and its numer- 
ical implementation. We present our numerical results 



for a stadium-shaped dot in section IIIIl We conclude in 
section IIVI Some technical details are provided in the 
appendices. 



II. MODEL AND METHOD 

Our setup consists of a rectangular sample of graphene 
of length L and width W :s> L, described by the Dirac 
Hamiltonian 

H = -ihva- ■ V + V{x) + C/(r), (1) 

where <t = {uxTUy) are Pauli matrices, v is the velocity 
of the Dirac quasiparticles, and U and V are gate po- 
tentials. The potential V{x) accounts for the metallic 
leads and is set to be V{x) -> — oo for x < and x > L, 
whereas V{x) = for < a; < L, so that the graphene 
sheet is tuned to the Dirac point in the region of interest; 
The potential U{r) defines the quantum dot. It takes a 
constant value U{r) = U inside the quantum dot, and 
C/(r) = outside the dot. 

In order to calculate the two-terminal conductance of 
the system, we apply Nazarov's matrix Green function 
formalismfSi which has been adapted to graphene by 
Titov et al^ Below, we summarize the essential elements 
of this method. Central object is a 4 x 4 matrix Green 
function G with matrix structure in the pseudospin space 
[corresponding to the Pauli matrices Ux and ay in the 
Dirac equation ([1])] and the "retarded-advanced" (RA) 
space. The matrix Green function is defined by the equa- 
tion 

(2) 

where e is the quasiparticle energy (set to zero in the re- 
mainder of the calculation) , and C is an additional param- 
eter. Following Ref. [13, we define the generating function 
J^((/)) as 

-^(^)-TrlnG-i|^^^^„(^/,^, (3) 

where the trace operation includes summation over RA 
and pseudospin indices, as well as integration over spatial 
coordinates. The generating function J- contains infor- 
mation about the full counting statistics of two-terminal 
transport through the systemi^ Here, we restrict our 
diskussion to the two-terminal conductance G, which is 
obtained from T through the equality 

^ = -250 S ' (4) 

where go = Ae'^/h is the conductance quantum in 
graphene. 

Denoting the matrix Green function for the graphene 
sheet without quantum dot {U = 0) as Gq, and writing 
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J"o = Tr In Gq ^ , we rewrite Eq. ^ as 

ST = T-Tq 

= Trln(l - C/Go). 



(5) 



The matrix Green function Go for the clean syste m ( with- 
out dot) has been calculated analytically in Ref. iSO. The 
explicit expression is rather lengthy and can be found in 
Appendix [X] Substituting the explicit expression for Gq 
into Equation (j4]) then gives 



G = 



9oW 
ttL 



SG, SG = -2go 



(6) 



4=0 



where the first term corresponds to the conductivity of 
clean graphene^'^'^ and the second term gives the cor- 
rection from the presence of the quantum dot, which is 
the focus of our calculation. 

The calculation of Ref. [s^ proceeds by expressing SJ- 
in terms of the T matrix of the quantum dot. Hereto, one 
introduces the Green function g for zero-energy quasipar- 
ticles in an infinite sample (i.e., with V = U = 0), 



5(r,r') 



as well as the function 



(T • (r ~ r') 



2iThv 



r-'|2 



G, 



eg 



Go - g, 



(7) 



(8) 



which, being a difference of two Green functions, is regu- 
lar if the spatial arguments coincide. With the standard 
Born series for the T matrixi^ 



one finds^ 



T = (1 - Ug)-'U, 
^^-Trln(l-rGrog), 



(9) 
(10) 



up to terms that do not depend on the counting field 4>. 
Equation pop is the basis for our further investigations. 

We will consider the limit that the size of the quantum 
dot is small in comparison to the length L. In that case, 
it is advantageous to expand Grog around the position ro 
of the center of the quantum dot. After a short algebraic 
manipulation, which is carried out in Appendix \^ one 
finds that one may replace the operators T and Greg in 
Eq. Uni) by matrices with elements T''" and G^^^, with 
/i, = 0, 1, 2, . . . and 



G^" 

^ rcg 



g^fj. grj./i^ 
1 



Greg(r,r') 



(11) 



r=r'=ro 



, , , d'^rd'^r'[{x-xo)-ia,{y-yo)f 



xT{r,r')[ix' - xo) + ia,{y' - ya)Y . (12) 



Note that each element G'^^^ is a 4 x 4-matrix with non- 
trivial operation on pseudospin and RA degrees of free- 
dom, whereas T^'^ acts in pseudospin space only and 
leaves the RA space untouched. With this replacement, 
the trace in Eq. (llOp is taken over the matrix indices ^ or 
V, the pseudospin degrees of freedom, and the RA degrees 
of freedom. 

Essentially, the matrix T^'' is the T matrix in a partial- 
wave expansion. In graphene, the partial wave expansion 
involves waves of angular momentum (7(/i + 1/2), with 
(7 = ±1 a pseudospin index and = 0, 1, 2, . . . the index 
of the matrix T'"^. The corresponding basis functions 
'0fc,/i,(T(r) are defined at a finite energy e = hvk only, 
where, for definiteness, we choose fc > 0. They are the 
solutions of the Dirac equation H^k^i^ai^) — e'4'k,fi,a{Y) 
in the absence of the potentials U and V . With the 
matrix notation '^k^l. = (V'fc,M,+i> V'fc,/i,-i), one has 



where J^, is a Bessel function and we used polar coordinates {r,9). Denoting the T matrix in the partial-wave basis 
by Tfj^^(k), one then finds that 

where we again refer to appendix |X] for details. In our discussion below, we will refer to the mode indices = 
0, 1, 2, . . . as s, p, d, . . . . 

It remains to describe a method to calculate the T-matrix T^^{k) for a specific quantum dot potential U{y). This 
is done through by first relating T to the scattering matrix S^^^^k) in the partial- wave basis, 

Sf^^{k) = Sf^i, - 2mT^_,y(k). (15) 

The scattering matrix Sf^y{k) relates the coefficients of incoming (— ) and outgoing (+) parts of the basis functions 
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(fT3| at energy e = hvk. Again using the shorthand notation ^'^^^ 



.;^), these are 




(16) 



where the i?;!^^ = ± ii^t are Hankel functions of 
the first (+) and second kind (— ), respectively. The 
scattering matrix S^^{k) is then defined through the 
asymptotic form of the solution of the Dirac equation 
H^pk{r) = /iffcV'fe(r) for r — > oo, which takes the general 
form 

Mr) - E + '^^\r)b,,) , (17) 

with, in 2 X 2 matrix notation, 

flfeM = ^Sfj^,y(k)bku- (18) 

To find Sfj,iy{k), we employ a variation of a method 
commonly applied to rectangular strips of disordered 
graphene.^** Hereto, we divide the graphene sheet in cir- 
cular slices flj-i < r < aj, j = l,...,N, with — 
and ajv so large that U{r) = for r > ajy- We then ob- 
tain the scattering matrix S{k) by solving two auxiliary 
scattering problems first. 

In the first auxiliary problem, we set the potential U to 
zero everywhere, except in the circular slice Oj-i < r < 
Uj. In this case, the wavefunctions can be expanded in 
incoming and outgoing partial waves for r < aj-i as well 
as for r > aj. The general solution of the Dirac equation 
in the regions r < aj-i and r > aj can be characterized 
by means of two reflection matrices pj and pj and two 
transmission matrices Tj and rj, such that p relates the 
coefficients of outgoing partial waves on the exterior to 
the coefficients of incoming partial waves on the exterior 
etc. If Sa = ttj — ttj-i is sufficiently small, these matrices 
can be calculated in the first-order Born approximation. 
Writing p^^ = Sp^^, p'^^ = r^^. = (5^^ (5t^^, and 

T^^ = (5^^ + hr'^y , this calculation gives 

= / ^^^^^ dr<)(r)t{/(r)*(;)(r), 

^P'^^ ^ 'Wvl dr<)(r)tt/(r)vl/(:)(r), 
^^'^^ = f dr<+)(r)t{/(r)*W(r), 

up to corrections of higher order in Sa. 

In the second auxiliary problem, we set U to zero for 
r > aj only. Defining Sj{k) to be the scattering matrix 
for this situation, we obtain the recursion relation 

S,{k) = p, + r;(l - S,^,{k)pr)-'S,^,{k)T,. (19) 



Together with the initial condition 5*0 = 1 and the equal- 
ity Sn = S this leads to the desired solution. 

Unitarity implies the relations Sp'j — ^^p]' ^'''j ~ 
—Stj, and Stj = —Stj^- Consistent with the Born ap- 
proximation, we may rewrite the recursion relation (jl9p 
as 

Sj{k) ^ Sj-i{k)[l + i6hj], (20) 

where 

iShj — S^j_i{k)6pj + Sp'jS.j-i{k) 

+ S]_,{k)6T;S,^^{k)+^T, (21) 

is a hermitian matrix. 

Although this concludes our formal description of the 
method, there are a few issues regarding the numerical 
implementation that we also need to discuss: 

(i) — For the final result, we have to take the low- 
k limit of the scattering matrix 5(fc), which, at first 
sight, may be problematic because the Hankel functions 
appearing in the Born approximation for the matrices 
dp, Sp' , St, and St' diverge in this limit. However, this 
problem can easily be circumvented by shifting the po- 
tential U{r) U{r) + hvk' for r < ajv, while at the 
same time increasing the quasiparticle energy according 
to k ^ k + k'. This means that we transfer contributions 
between the free propagation and the perturbation. We 
then solve the scattering problem inside the disk starting 
from "free" Dirac fermions with energy hv{k + k'), which 
remains finite in the limit fc — >■ 0. The wavefunction out- 
side the disk, where free electrons have momentum fc, 
can be found by matching to the wave function inside 
the disk, and contains the information of the scattering 
matrix we need. 

For the low-fc-expansion of the scattering matrix, we 
then find an expression of the form 

- + 5;,(fc')fc^+''+' + 0(fc''+^+'), (22) 

where S'{k') depends on ojv, k' , and the scattering ma- 
trix S^ik') obtained for the modified problem where we 
have replaced fc by fc' -I- fc and safely can take the limit 
fc ^ 0. A detailed derivation is given in Appendix |B] 
Note, that the expansion in Eq. (1^^ gives us precisely 
that order in k, that we need in order to obtain the ma- 
trix T'"^ in Eq. ([Ti)) . 

The wavenumber k' of the quasiparticles can be cho- 
sen arbitrarily, and we have verified that the result of our 
calculations do not depend on the choice of k' . However, 
one may exploit this freedom for a wise choose of k' . It 
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should be not too small, because otherwise the Hankel 
functions in the spherical incoming and outgoing waves 
become too large, indicating the fact that particles are 
repelled from the origin and therefore the Born approxi- 
mation looses its applicability. On the other hand, if the 
wave vector k' is chosen too large, the effective potential 
U + hk'v gets large, too, and the thickness da of the slices 
has to be decreased. 

(ii) — It is sufficient to choose the radius gq of the 
first slice such, that U{r) is uniform for r < oq. For this 
situation an exact solution is available. 

(iii) — In order to guarantee numerical stability during 
multiplication of scattering matrices, we change Eq. ((20)) 
to 

which is valid up to corrections beyond the accuracy of 
Born approximation, but ensures that the scattering ma- 
trices remain unitary at all times. 

(iv) — In the calculation of the scattering matrix, we 
keep a total of M modes. After the calculation of the full 
scattering matrix S, we then further truncate the scat- 
tering matrix upon calculation of the determinant in Eq. 
pO)) . keeping matrix elements S^^ with /i, i/ = 0, 1, M' 
only. Keeping a large number of modes M in the calcu- 
lation of the scattering matrix is important in order to 
properly resolve the dynamics inside the quantum dot. 
The number of modes M' required for the calculation of 
the conductance depends on the ratio of the size L of the 
graphene sheet and the size of the quantum dot, and can 
be kept small if this ratio is large. The numerical calcu- 
lation of the conductance has converged, when the result 
does no longer change upon increasing M, M' or N. For 
the results of the conductance shown here, we kept chan- 
nels up to d-wave {M' = 2) only, which we found to be 
sufficient for the parameters chosen in our calculation. 



III. TWO-TERMINAL CONDUCTANCE WITH 
STADIUM-SHAPED QUANTUM DOT 

A. Two-terminal conductance 

As a nontrivial application of the formalism laid out in 
the previous Section, we have calculated the two-terminal 
conductance of a stadium-shaped quantum dot in an oth- 
erwise undoped graphene steet. The quantum dot is 
placed halfway between the contacts. It is character- 
ized by the radius R of the circular pieces and the length 
2a of the linear segment, see Fig. [TJ The gate voltage 
U{r) = [/ for r inside the quantum dot, and U{r) = 
otherwise. In our calculations, we have included contri- 
butions up to d-wave in the T matrix, whereas as many 
modes were used for the calculation of the scattering ma- 
trix S as were necessary to reach convergence. (See the 
discussion at the end of the previous section.) 
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FIG. 2: (Color online) Conductance versus gate voltage U 
for an undoped graphene sheet with a gate-defined stadium- 
shaped quantum dot. The parameters of the numerical cal- 
culation are 2a/R = VS and L/R = 5. The curves corre- 
spond to the numerical calculation using the method of Sec. 
im with contributions to the T matrix up to d-wave (blue, 
solid), and to a direct numerical solution of the Dirac equa- 
tion (red, points; data taken from Ref . |20|) . 



In Fig. [21 we show the two-terminal conductance as a 
function of the gate voltage U. The conductance exhibits 
several resonances due to the formation of quasi-bound 
states in the stadium. The parameters for the data in this 
figure are chosen to be the same as in Ref. [13, so that 
the results can be compared. The excellent agreement 
between both curves concerns the position as well as the 
shape of the resonances. We attribute small deviations 
in the resonance heights to corrections from higher angu- 
lar momentum that were not taken into account in the 
calculations. Such corrections disappear upon increasing 
the ratio of system size vs. dot size L/R (see discussion 
below). 

The effect of inclusion of successive angular momen- 
tum channels in the T matrix (while keeping essentially 
all angular momentum channels in the calculation of the 
scattering matrix 5*) is illustrated in Fig. [31 where, in 
the upper panel, we show a close-up of the data of Fig. 
[2] with one, two, and three angular momentum channels 
included in the final step of the calculation (i.e. with 
truncation at M' = 0,1,2). This figure clearly demon- 
strates that each resonance has contributions from more 
than one angular mode, as is to be expected for a chaotic 
quantum dot. In particular, we find that all resonances 
have a non- vanishing s-wave contribution, so that the res- 
onance position can be extracted from the s-wave channel 
solely. For the data shown in Fig.[2l one needs to include 
the p-wave contribution in order to obtain the correct 
resonance shape. The effect of the d-wave contribution 
is small, which we find remarkable, because the ratio of 
dot size 2(i? -|- a) versus contact size L is « 0.75, which 
is not small in comparison to unity. The lower panel of 
Fig. [31 shows conductance data for M' = 0, 1, 2, showing 
that the convergence with respect to M' quickly improves 
upon increasing L/R. 

We conclude, that most information about the reso- 
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FIG. 3: (Color online) First four conductance resonances: 
The conductance has been calculated using Eq. (|10p with 
the T matrix truncated after the s-wave (blue, dashed), p- 
wave (red, dotted), and d-wave angular-momentum channel 
(green, solid), corresponding to M' = 0, 1,2, respectively. In 
the upper panel, we set L/R = 5, while in the lower panel, 
L/R= 10. The ratio a/R remains the same as in Fig. [2] 



nances that is relevant for transport at large or mod- 
erately large L/R is encoded in the contributions from 
small angular momentum, and is therefore stored in a 
small set of parameters. That fact is a strong indication 
of the power of the method employed here, in comparison 
to a direct simulation of the Dirac equation. 



B. Lineshape 

The main advantage of the present method over the 
direct numerical solution of the Dirac equation is that 
it allows one to extend the conductance calculations to 
the regime L ^ i?: The complete dependence on the 
length L is encoded on the Green function Greg, which 
is known analytically. It is only in the regime L ^ R 
that resonances are well separated and can be charac- 
terized individually. The calculations of Ref. [13, on the 
other hand, were limited to the regime L < 5i?, where 
resonances were still strongly overlapping. (Reference [20l 
also considered the case L sa 8R for W = L. However, in 
that case the conductance is strongly influenced by the 
finite width W of the graphene sheet.) 

In Fig. 2] we show how the first and second resonance 
behaves upon increasing the contact size L at fixed dot 
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FIG. 4: (Color online) Dependence of the first and second 
resonance on the ratio of dot size versus system size R/L. 

size R. The plots illustrate that the resonance width 
shrinks, while the height saturates. The asymmetry of 
the line shape disappears in the limit L/R — > oo, consis- 
tent with the expectation that s-wave scattering domi- 
nates if L ^ i?. Similar behaviour is found for the other 
resonances. 

In order to quantitatively analyze the lineshapes, we 
note that each resonance is characterized by a diver- 
gence of the T matrix. With a suitable parameterization 
of these divergencies an explicit expression for the reso- 
nance lineshape can be obtained. Hereto we introduce 
the dimensionless variable e = UR/hv. Then, close to a 
divergence at eo , the divergent part of T^'' is of the form 

Rf^" 

T^^^ie) ~ (24) 

e — eo 

where the matrix R^^'^ contains the information about the 
resonance shape. The matrix elements R^'^ are related by 
hermiticity, i?^"^ ~ (7?"^)^ and time-reversal symmetry, 
i?^'' = ay{R'^^^Y^ay. Moreover, for the specific problem 
we consider here, inversion symmetry and reflection sym- 
metry impose further constraints on R, which allow us 
to parameterize the matrix as 

-ilspUx IpO-o ilpd<7x , (25) 

lsd<70 —ilpdC^x IdC^O / 

where ctq is the 2x2 unit matrix. Due to Kramers degen- 
eracy, R can be decomposed into two identical submatri- 
ces, which, for a generic resonance, are of unit rank. (In 
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a chaotic quantum dot level repulsion ensures that this 
condition is always fulfilled.) This gives the additional 
constraints 
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/2 7 7 72 7 7 



(26) 



on the parameters in Eq. (I25t . Note, that , Ip, and 
Id have the dimension of a length, length^ , and length'"^ , 
respectively. 

Substituting Eq. ([25]) into Eq. pU)) and making use 
of the condition (PB)) . we find that the line shapes are 
described by the Fano resonance formula— 



|2(£-eo) + grr 
■ 4(e-eo)2+r2 



(27) 



Here, Gnr is the non-resonant conductance, T is the res- 
onance width, eo is the resonance position, and q is the 
complex "Fano parameter" . After subtracting the back- 
ground conductance Gnr, we rewrite Eq. (I?71) as 



(5G = go 



(i + Aa{e - ep) 
4(e-eo)2 + r2' 



(28) 



where /3 = {Gnr/9o){\q? - 1)^^ and a = (G„r/go)Reg. 
For a = the resonance has a Breit-Wigner shape; Non- 
zero a is responsible for an asymmetry in the lineshape. 
When we express the resonance parameters through the 
entries of R, we also use asp = sign(Zsp) etc., which is not 
fixed by Eq. (HH): 



eo = 

r = 

13 = 







£0 + 


24L2 


Is 


9 7 

TT /„ TT 


2L 


+ 8L3 + 



960L4 




(29) 
(30) 



180L6 
25iyy|crgrf7r'^ 



18L6 
13549/^7r6 
28800L10 

'Ipld^pd 



90i8 



(31) 
(32) 



Note, that the shift ~ of the resonance position, 
as well as the asymmetry a is absent, if the resonance 
consists of a single angular momentum component only. 

Now consider the situation, where we fix the product 
UR of the gate voltage quantum dot and the quantum 
dot size — i.e., we look at a fixed resonance — and in- 
crease the length L of the graphene sheet. Then, as long 
as Is 0, at sufficiently large L, the conductance will be 
determined by its s-wave contribution, so that the reso- 
nance line shape reduces to the Breit-Wigner form. In 
this limit, the height for SG approaches the constant uni- 
versal value Sgo/ir'^, and the width scales as F = Is/^L, 
as was found previously by Titov et aZ— 



C. Comparison Disc - Stadium 

We now summarize the main qualitative difference in 
the conductance resonances for stadium-shaped quantum 
dots, which were considered numerically here and in Ref. 
|20| . and disk-shaped dots, which were considered ana- 
lytically in Ref. HO and numerically in Ref. llOj. The 
main difference, which is illustrated in Fig. |31 is that 
generic resonances for the stadium dot have contribu- 
tions from more than one angular momentum channel. 
This means that the resonance lineshape changes from a 
Fano shape for moderate ratios L/R to a Breit-Wigner 
lineshape for large L/R, as seen in Fig. 01 In the limit 
L/R ^ oo of isolated resonances, the height of the con- 
ductance resonances takes the universal value Sgo/n^, 
whereas the width of all resonances scales proportional to 
R/L. Though we considered a stadium-shaped dot only, 
we believe, that the features found here are generic and 
shared by all chaotic dots. 

This is to be contrasted to the situation of a circular 
quantum dot, where different angular momentum chan- 
nels do not mix, and the resonances are of pure type 
(s-wave, p-wave,. . . ). The resonance lineshape is always 
of Breit-Wigner form and approaches a constant height 
upon taking the limit L/R — > oo. The asymptotic reso- 
nance height depends on the angular momentum channel, 
with a height Sgo/ir'^ for s-wave resonances {I = 0) and 
< 2gQ for higher angular-momentum channels {I > 1). 
The resonance width depends on the angular momentum 
channel, F cx {R/L)^'+\ 

We close this comparison with a comment on the re- 
sults of Ref. [l^. There it was found that the resonance 
height of the chaotic dot goes to zero upon taking the 
limit of large L/R. This observation referred to a setup 
with aspect ratio W/ L = 1 and periodic boundary con- 
ditions along the transverse direction. For aspect ratio 
W L, only the mode with zero transverse momen- 
tum substantially contributes to transport.— Therefore, 
s-wave scattering does not affect the conductance, re- 
sulting the observed suppressed height of s-wave conduc- 
tance resonances. We verified that the same phenomenon 
occurs in our calculations, by evaluating the regularized 
Green function Greg for finite width. This explains, why 
the resonances in Ref. i2^ completely disappear in the 
limit of large L/R, while in our investigation for aspect 
ratio W ^ L the height of the resonances remains finite. 



IV. CONCLUSION 

In this article, we investigated the resonances of the 
conductance of a graphene sheet with a chaotic quan- 
tum dot. Using a numerical implementation of the ma- 
trix Green function method of Ref. [s^, we were able to 
study the behaviour of the resonances in the limit of well- 
separated resonances. This essential limit could not be 
reached in the original treatment of the problem. 

As was proposed in Ref. (2^ the resonances of the 
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chaotic dot behave significantly different compared to the 
case when the gated region is circular. While the circu- 
lar geometry does not allow for mode mixing, so that 
all resonances are Breit-Wigner resonances with a well- 
defined angular momentum, in the stadium dot, due to 
its non-integrable dynamics, all resonances have contri- 
butions from all scattering channels. The presence of 
mixed angular- momentum modes is responsible for an 
asymmetry of the line shape, described by the Fano res- 
onance formula. In the limit of very well separated res- 
onances, corresponding to the limit in which the size R 
of the quantum dot is much smaller than the distance 
L between the metallic contacts, for the chaotic quan- 
tum dot all resonances are dominated by the lowest (s- 



wave) angular momentum component, and recover the 
Breit-Wigner form. In contrast, for the disk-shaped dot, 
parametrically narrower resonances for higher angular 
momentum channels (p-wave, d-wave, etc.) persist in 
the limit of large L/R. This establishes clear signature 
that distinguishes regular and chaotic dynamics of gate- 
defined quantum dots from the conductance resonances. 

We gratefully acknowledge discussions with J. H. Bar- 
darson and M. Titov. This work is supported by the 
Alexander von Humboldt Foundation in the framework of 
the Alexander von Humboldt Professorship, endowed by 
the Federal Ministry of Education and Research (PWB) 
and by the German Research Foundation (DFG) in the 
framework of the Priority Program 1459 "Graphene" . 



I 

Appendix A: Matrix Green Function 



In this appendix, we give explicit exp ressions for some of the matrix Green function appearing in Sec.|TIl We follow 
the supplementary material of Ref. l30l . but our results for the d-wave channel go beyond that reference. 

Using Pauli matrices t^, Ty, and for the retarded-advanced (RA) degree of freedom, and with tq for the 2x2 
unit matrix in the RA grading, the explicit expression for the matrix Green function Gq is 

Go{x,x';y) = -^V{x)AArAaAV~\x') (Al) 















2 




n 



with A= „ „ ] <S)(Tz+ { n 1 ^ctq and 



r _ f icosh{(j)y/2L) smh{^y/2L) \ 

- sinh((/.?//2L) -zcosh(0y/2L) )^^°' ^^^> 

Act = Tfj® \ '^'"l'f(^+=^j'+*a)/2il sin[iT{x-x^+iy) /2L\ \ ^ {Ai) 
\ s\n\'n{x~x' ~iy) /2L sin[7r(2:+a:' —iy)/2L J 

y{^) = % . % )®^o. (A4) 

\ I COS jj- z sm 1^ / 

Ths s-wave contribution of the regularized Green function can be found as 

GtL{xf)) = 1™ [Go{xo, x'\ y) - g{xo - x', y)] 

Contributions from higher angular momentum modes are obtained by keeping higher terms of the Taylor series. 
Therefore, we write the regularized Green function as 

aeg(r, r') = ^^[(^ - ro)^]''G,eg(ro, ro)[^'(r' - ro)]"- (A6) 

Here, ^ acts to the right, on the first argument of the Green function, and V' acts to the left, on the second argument 
of the Green function. We can simplify this expression by using the equations of motion of the regularized Green 
function, 

-ia'^Greg^O, Greg(-«V'cr) = (A7) 

from which immediately follows, that 

5j/Grog(r,r') = iazdxG^c^{v,v'), 5j,'Grog(r, r') = dx'Gxc^iiy , ^(yz)■ (A8) 
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This enables us to rewrite Eq. (jA6l) as 
with 



G,cg(r,r') = ^[(^ - 2^0) + ^^.{y - V^WG^^M^' - xo) - ia.{y' - yo)]^ (A9) 



G^:^^d:id:,G,,^{ro,ro). (AlO) 

We now give the exphcit expressions for G(4g d-wave order. In order to keep the expressions short, we focus 

on the case, where the dot is placed in the middle of the sample {xq — L/2), 

/^sd/c 
reg 





f^pdj dp 

The element G^^g is given by Eq. (|A5[) . Since V{xq) has no matrix structure in pseudospin space, it commutes with 
the T-matrix, and therefore does not play a role for the generating function. 
Upon inserting (IA9|) into the generating function (Eq. (|10l) ). we find 

(5J" = lndet[l-TGrcg], (All) 

where Grog is a (infinite) matrix with entries G(^(,g {^,v — 0, 1, ...) and the matrix T contains the elements T^^^ 

T^^- j £vSv'j,[{x-x^)-ia,{y-y^)rT{v,v')^[{x'-x,,)+ia,{y'-y^)r. (A12) 

This is Eq. p2|) of the main text. In the limit that the size R of the quantum dot is much smaller than L, this expansion 
is convergent, and it is a good approximation to limit the number of angular momentum "channels" involved in the 
expansion. In the simplest case, one takes into account the s-wave contribution fi = v = only. This limit was 
considered in Ref. [st^ The accuracy can be improved by including contributions from higher angular momentum. 
The expressions derived above allow one to go up to the d-wave contribution. 

In order to relate the object T^"^ of Eqs. (fT^ or (jA12[) to the T matrix in the partial- wave expansion, it is instructive 
to compare the plane- wave and circular- wave basis sets of eigenstates of the free Dirac Hamiltonian Hq = — icrV. The 
plane wave basis consists of eigenstates, which are labeled by their wavevector k = (k cos 9k, ksinOk), 

Mr)-"^[X)- (A13) 

Here and below, we restrict ourselves to positive energy solutions (conduction band) only. In the implementation, this 
can be achieved by choosing the gate potential of the dot to be negative. The circular- wave states are combined eigen- 
states of the Hamiltonian and the total angular momentum Lz + ■^o'z, and are therefore labeled by their wavenumber 
k and their half-integer angular momentum quantum number m, 

^,„.„.;^,-f y (A14) 

V 47r \^ e'2zsgn(TO)J|„^^i|(fcr) y 

Here, Ji, is the Bessel function of z^-th order. The 2x2 matrix ^f/j^ of the main text is related to the ipkm as 
^fe/i = (V'fc,/i+i/2: *'0fc,-/i-i/2)- The basis change between plane waves and circular waves is expressed through the 
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equation 

V'k(r) = ^'^-^'e-^^^-^^'^V^fonW. (A15) 

m 

For the scattering problem in the infinite graphene sheet, we place the center of the quantum dot at the origin. We 
may express the T-matrix in the plane-wave basis through the T-matrix in real space as 

T(k,k') ^^J (frd^r'{l e-''^'')T{r,r') ^ Jg,^ ^ exp[ifc(x'cos6lfc + y'sin^;, - xcos6lfe - ?/sin6lfe)]. (A16) 

On the other hand, the same matrix elements can be also written in spherical wave basis, 

T(k,k') = ^ V(-z)l"-5lil™-^le*(»-3)«^e-*(™-^)''^r„™(yfc) (A17) 

nm 

The partial- wave T matrices Tmn introduced here and T^^ of the main text are related as 

rp ^ ( 7)i+l/2,i/+l/2 i7)i+l/2,-i/-l/2 j (A18) 

Comparing Eqs. (jA16l) and (|A17I) . we obtain a relation between the T-matrix in spherical wave basis and the T-matrix 
in real space. For example, for the s-wave channel we find 



2tt f Ti i(k) iTi _i(k) 

I 2-2 2'2 

-iT_i i(k) T_i _i(fc) 



= i y d^rd\'T{r,r') + Oik), (A19) 



nn , 47r / Ti i(k) iTi _i(k) 



so that 

fc^o k \ -iT_i i{k) T_i _i{k) 

\ 2'2 2'2 

= lim^Too(fc), (A20) 
fc->o k 

where, in the last line, we used the 2x2 matrix notation with integer indices used in the main text. The remaining 
identities in Eq. (fM)) follow similarly. 



Appendix B: Low-fc-Limit of Scattering matrix 

In this appendix, we present the details how to determine the low-fc limit of the scattering matrix. As explained 
in the main text, we modify the scattering problem, such that we solve the scattering problem for particles with 
wavenumber ki = k + k' in the potential U' = U + hk'v. We then need to relate the scattering matrix S' for this 
problem, in which U' is set to zero outside a certain cut-off radius a to the scattering matrix S for the original problem, 
in which U, not U' is set to zero for r > a. 

In order to relate these two scattering matrices, we need to calculate the scattering matrix of a potential step, in 
which the potential is zero for r > a and equal to it = hk'v for r < a. Hereto, we expand the wavefunctions for r < a 
and r > a in terms of incoming and outgoing circular waves with wavenumber k' + k and fc, respectively, 

^a-ir) = ^(L+V'fc+fe.,™^+(r) + L;;^Vfc+fe/^„,_(r)), 

m 

^a+(r) = ^(i?+V'fc,™,+ (r) + i?^Vfc,™,-(r)), (Bl) 

m 

where the scattering states V'/c,m,± are obtained from the basis states ipkm in Eq. (IA14[) by the replacement J{kr) — > 
H^{kr). The coefficients i?^ and are related through the transfer matrix T, 

( ) " k + k',a)(^ ^ , (B2) 
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which is easily calculated from continuity of the wave function at r = a, 



T„,{k,k + k',a) = — fmik,k + k',a)^k{k + k'), (B3) 

where the 2x2 matrix Tm{k, k + k' , a) has entries 

fm(k, k + fc', a)n = H\-^^^,{ka)H\+^_, (ka + k'a) - h\^^_, {ka)H\^^^, {ka + k'a), 
r™(fc, k + fc', a)i2 = Hl^^^^,{ka)HlJ^_, (ka + k'a) - ^ (fca)irf^|^ ^ {ka + k'a). 

The other entries are given by the relations 7m(fc, k + fc', 0)21 — —%n{k, k + k' , a)|2, Tm(k, k + fc', 00)22 = —%n{k, k + 
k',a)ii. The transfer matrix is converted into reflection/transmission coefficients using the standard relations 

P = Tl2(7^2)"\ 

T - (r22)-\ 

t' = Til - Tl2(722)~^721, 
P' = -iT22)-'T2l. 



(B4) 



Upon taking the limit fc — >■ 0, we find 



■ „2|m| 

fcl™l +C'(fcl"l+i), 



™ (|n^|-i)!2l"l~l Hl-l^{k'a) 
Hl+1_,{k'a) 

We can now relate the scattering matrices S' for particles with energy k' in the potential U' = U + hk'v and ?7' = 
for r > a to the scattering matrix S for particles with energy k in the potential U, with U — for r > a, 

S' = /9 + r'(l-5y)-iS"r. (B6) 

Inserting the above expansions, we get an expression of the form 

Snm{k) = 5nm + S^^l A^l^l + I"! + O (fcl"l + l™l + l ) , (B7) 

coinciding with Eq. (|22p in the main text. The coefficient S*^^-* depends on a, k' and S' . 
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